Supercurrent and dynamical instability of spin-orbit coupled ultracold Bose gases 
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We investigate the stability of supercurrents in a Bose-Einstein condensate with ID spin-orbit 
and Raman couplings. The consequence of the lack of Galilean invariance is explicitly discussed. 
We show that in the plane wave phase, characterized by uniform density, the supercurrent state 
can become dynamically unstable, the instability being associated with the occurrence of a complex 
sound velocity, in a region where the effective mass is negative. We also discuss the emergence of 
energetic instability in these supercurrent states. We argue that both the dynamical and energetic 
instabilities in these systems can be generated experimentally through the excitation of the collective 
dipole oscillation. 
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A characteristic feature of superfluids is their capa- 
bility to support a current flow (supercurrent) without 
dissipation [HQ- In a Galiean invariant, uniform config- 
uration the supercurrent does not decay if the velocity of 
the fluid is smaller than the critical velocity fixed by the 
famous Landau criterion 0, Q v CT = mm. e(p)/p where 
e(p) is the excitation energy of the elementary excitation 
of the system carrying momentum p. Landau's instability 
has an energetic nature, being associated with a negative 
value of the excitation energy. Non-uniform supcrfluid 
systems, like ultracold Bose-Einstein condensates in op- 
tical lattices, are known to exhibit also dynamic instabil- 
ities. Both the energetic and dynamic instabilities have 
been the subject of intense theoretical and experimental 
investigations in ultra cold atomic gases |5l4l2|. 

The recent realization of spin-orbit coupled gases 
[l7T ] is opening new perspectives in the study of superfluid 
phenomena. These systems lack Galilean invariance and 
show its consequences also in uniform density configu- 
rations. In particular it has been recently shown that 
in these systems the usual Landau criterion for stabil- 
ity, which applies to the motion of an impurity in the 
medium, cannot be used to determine the stability of 
configurations carrying a supercurrent, the correspond- 
ing critical velocities being dramatically different 3, 1|| . 

In the present work, we discuss the dynamical as well 
as the energetic instabilities exhibited by supercurrent 
states in an ultracold Bose-Einsetin condensed gas with 
a type of spin-orbit coupling already realized in experi- 
ments [13, llTf • We first discuss the connection between 
the lack of Galilean invariance and the fact that the cur- 
rent is not conserved in such systems. We then determine 
the supercurrent state by adding a Lagrange constraint 
to the Hamiltonian and explore the instability conditions 
by calculating the corresponding frequencies of the ele- 
mentary excitations. The critical values of the velocity 
are calculated as a function of the Raman coupling SI. 
Finally we show that the collective dipole oscillation, in 
the presence of a harmonic trap, is well suited to investi- 
gate the effects of instability, and make first comparisons 
with a recent experiment [17] . 



Hamiltonian and lack of Galilean invariance. — We 
consider a two-component Bose gas with spin-orbit cou- 
pling whose Hamiltonian is given by dUEll 
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where Vv,p is an annihilation operator of a particle with 
spin a and momentum p. The Hamiltonian ([T]) can be de- 
rived by applying a local (ir-dependent) rotation in spin 
space to the Hamiltonian defined in the laboratory frame 
where two detuned ^spin-polarized laser fields are coupled 
to the system 



13 



225. The spin-orbit coupling strength 
fco is fixed by the momentum transfer of the two lasers, 
while the Raman coupling Q is determined by the laser 
intensity. The matrices Oi with i = x,y, z are two-by-two 
Pauli matrices. The single-particle part of the Hamil- 
tonian (JTJ) has two branches, the lower energy branch 
exhibiting two minima if fl < ^k^jlm. We assume an s- 
wave interaction with a common intra-species interaction 

(.9tt = 9U = 9)- 

"Hint =^ X! E V^p^paVv^Vv.p! 

ff=td Pl+P2=P3 + P4 

+ 9 ~f E ^pXps^p^T,^ (2) 

Pl+P2=P3+P4 

where V is the volume of the system. Despite the pres- 
ence of the laser fields in the laboratory frame, the Hamil- 
tonian ([T]) in the spin rotated frame is translationally 
invariant, the commutation relation [H, P] = being 

exactly satisfied where P = £ p p (vjpV'tP + ^Ip^-lp 

The Hamiltonian (|T|) lacks 



is the momentum operator, 
however Galilean invariance 
the fact that the current operator 



This follows from 
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satisfying the equation of continuity does not coincide 
with the momentum, due to the presence of the spin term. 
Furthermore, due to the presence of the Raman coupling 
in ([1]), the current operator does not commute with H. 

Supercurrent state. A natural way to construct a 
supercurrent state is through the determination of the 
ground state of the Hamiltonian 



H v = H-v-P 



(4) 



obtained by adding to H a Lagrange multiplier term pro- 
portional to the total momentum P. Using the operator 
P rather than the current operator J in ((4]) is crucial in 
order to ensure that the ground state of H v is also an 
eigenstate of the original Hamiltonian (pj. This is guar- 
anteed by the commutativity of P with %. 

In the present work, we focus on the plane- wave phase, 
in which only a single momentum state is macroscopically 
occupied and the density of the system is uniform. In the 
plane wave phase the macroscopic wavefunction of the 
condensate takes the form 
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where n is the density of particles, and and ki are 
variational parameters to be determined by minimizing 
the energy with respect to H v . The angle determines 
the spin struture of the wave function. In particular the 
average spin polarization of the gas is given by (a z ) = 
cos 2$, while ki determines the value the momentum P, 
equal to iVki, with N the number of particles. 

The energy per particle E v = (H v )/N, calculated on 
the state ip v , is given by 
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Setting the derivative of E v with respect to ki equal 
to 0, one obtains the relationship 



k : = (ko cos 20,0,0) 



(7) 



in which case we have 
cos 20 = 



y/kf + (mft/2fc ) 2 ' 



(8) 



independent of the value of the interaction coupling con- 
stant. In © we have chosen sin 20 > and the sign of 
cos 20 to be the same as that of k% , which corresponds to 
choosing the lower branch configuration of the system. 
From (0 and (jHJ, one can express the velocity of the 
supercurrent state as a function of k\ as 
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The stationarity conditions for k\ and give, in general, 
three solutions (see Fig. QJc) andQJd)), one correspond- 
ing to the lowest energy state (ground state of H v ) where 
v and k\ have the same sign, a second metastable solu- 
tion corresponding to a local minimum where v and k\ 
have opposite sign, and finally a solution corresponding 
to a local maximum of E v . Below, we discuss the stabil- 
ity of these supercurrent configurations. To this purpose, 
in the following we focus on the case g = g^, which is 
a good approximation in the available experiments with 
Rubidium atoms 13 , 13] • It is straightforward to extend 
the formalism to the case g ^ g^. 

Dynamical Instability. — The dynamical instability is 
characterized by the occurrence of positive imaginary 
components in the excitation spectrum. We find the con- 
ditions for the dynamical instability by looking at the 
dispersion of the elementary excitations and in particu- 
lar at the sound velocity. The dispersion of the elemen- 
tary excitations is obtained by investigating the poles of 
the Green's function within the Bogoliubov approxima- 
tion [23|], taking the condensate wavefunction ([5]) with 
the value of k\ and determined by the solution of 
and (JSj> . (The dispersion can be equiyalently obtained us- 
ing the hydrodynamic formalism [22j.) In the long wave- 
length (sound) limit we can write the dispersion in the 
linear form uj w c±|g|, where q is the momentum shift rel- 
ative to k\ and c+ (c_ ) is the sound velocity when q > 
(q < 0). After some striaghtforward algebra we obtain 
the result 



between the value of ki , the angle and the velocity v. 
Recalling that cos 20 gives the average polarization one 
immediately finds that the total current of the state ip v is 
given by the expected relationship (J) = Nv, revealing 
that the velocity of the flow in the state ip v coincides 
with the velocity v entering the Lagrange constrained 
Hamiltonian % v . In the following, we assume that v, 
and hence ki, be oriented along the x-direction, and to 
simplify the notation we set k\ >x = k\. 

Similarly, setting the derivative with respect to equal 
to 0, one obtains an equation for 0. This equation takes 
a simple form when the interaction is isotropic (g = g^i) , 



where 



k\ 
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is the lower branch eigen-energy of the single-particle 
part of the Hamiltonian (fTJ. The sound velocity (fT0|) 
develops an imaginary part when the inverse effective 
mass d 2 t-jdk\ becomes negative 24|. The region of 
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dynamical instability is marked as red dotted lines in 
Fig. [T] (a) and (b). A similar mechanism of dynamic in- 
stability is exhibited by quantum gases in the presence of 
periodic potentials [5l4l2| where the negativity of the ef- 
fective mass is caused by the band structure of the single- 
particle excitation spectrum Q . There is however a deep 
difference between the two cases. In fact, differently from 
the case of gases in a periodic potential, the density of 
our systems is uniform and the violation of Galilean in- 
variance is caused by the presence of the spin term in the 
current operator ([3]). 

The spin structure of the order parameter of the sta- 
tionary solution of T~i v coincides with the one of the non- 
interacting model. In particular, the velocity ((9]) of the 
supercurrent state is identical to the velocity de-/dk\ of 
the single-particle state with momentum k\ . This corre- 
spondence enables us to understand the stationary solu- 
tions of E v in terms of the single-particle eigenstates. As 
shown in Figure Q] (a) and (b) , there are up to 3 single- 
particle eigenstates for a given velocity v. In Figure [2(c) 
to (f), we plot E v with the condition © for the spin 
structure of the wavefunction, as a function of k\ for dif- 
ferent values of v, as indicated in Figure Q] (b) . When v 
is positive, the energy minimum state with positive k\ is 
the lowest energy state, while the energy minimum state 
with negative k\ (if exists) is a mctastablc state. One can 
notice that the local minimum energy states of H v cor- 
respond to single-particle eigenstates having eigenenergy 
e-(ki) with positive curvature. Conversely, local max- 
imum energy states of T-L v correspond to single-particle 
eigenstates having eigenenergies with negative curvature, 
and are dynamically unstable. 

Energetic Instability. — The energetic instability oc- 
curs when the system develops a negative excitation en- 
ergy. As before, we focus ourselves to the case where v 
is oriented along the x direction and we consider equilib- 
rium configurations with positive k±. When v is small, 
both the lowest energy state and the metastable state 
of local minimum energy have excitation spectra every- 
where positive. However, above a certain critical value of 
v, the excitation spectrum starts taking negative values. 
Figure [2] plots the critical velocity as a function of O. 
For the interaction parameter^we take gn = 0.48/cQ/2m, 
which is the value used in [22fl and is relevant in the 
presently available experiments [17[ . The blue solid line is 
the critical velocity calculated assuming that the system 
is in the minimum energy state (v and k\ have the same 
sign), which has been already calculated in [19]. The 
green dashed line instead corresponds to the metastable 
configuration (v and k\ have the opposite sign). The 
metastable state does not exist below the red dotted line, 
corresponding to the onset of the dynamical instability, 
which approaches the value —ko/m as Q — > 0. 

We find that, when is small, the instability of the 
lowest energy state and the metastable state is due to the 
roton part of the excitation spectrum. At large f2, the 




FIG. 1: (a) Single-particle eigen-energy and (b) correspond- 
ing speed as a function of k\ when Vl — 1.7kg /2m. (c) - (f) 
E v for four different values of v (indicated in (b)), as a func- 
tion of fci . Local minima of E v are marked by orange points 
and local maxima by black points. In (a) and (b), the dy- 
namically unstable region is shown in a red dotted line, the 
region which corresponds to the metastable states in green, 
while the lowest energy states of E v are in blue. The collec- 
tive dipole oscillation occurs around one of the minima of the 
single-particle spectrum, as described in (a). 




FIG. 2: Critical velocity for the energetic instability of the 
supercurrent, in units of ko/m, plotted against Q, in units 
of fcg/2m, when the condensate has a positive momentum. 
The blue solid line corresponds to the lowest energy state, 
and the green dashed line to the metastable state, when 
gn — 0.48fco/2m. The red dotted line provides the onset 
of dynamical instability. 
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instability is instead caused by the phonon branch. The 
transition between the roton instability and the phonon 
instability for the lowest energy state occurs around 
S7 ps 4.6fcg/2m, where one can see a kink in the fig- 
ure. In the absence of the supercurrent, the phase transi- 
tion between plane- wave phase with nonzero momentum 
and the condensate with zero momentum takes place at 
SI = 4.0fcQ/2m, and when the condensate occupies the 
zero momentum single particle state, the excitation spec- 
trum does not exhibit any roton structure. Nevertheless 
in the interval 4.0k 2 /2m < SI < 4.6k 2 /2m the instability 
of the supercurrent occurs due to the emergence of ro- 
ton excitations, which is characteristic of the plane-wave 
phase. The transition between the roton and phonon 
instabilities for the metastable state occurs at smaller 
values of the Raman coupling, around S7 ps 2.0k 2 ) /2m. 
The maximum (in magnitude) velocity compatible with 
the existence of a metastable state (the red dotted line in 
Fig. [2j also see Fig. Q] (e)) approaches zero as S7 reaches 
SI = 4.0/c§/2m from below; above S7 = 4.0fcg/2m there is 
no metastable state. If we introduce a small difference 
in the interaction parameters g ^ g^, the overall struc- 
ture does not change in a significant way, except at small 
values of £7 where the critical velocities vanish at a finite 
value of SI, corresponding to the transition to the stripe 
phase 0. 

One can also discuss the instability in terms of the 
momentum kick (or momentum excess) that one can pro- 
vide to the condensate, initially at equilibrium. In Fig- 
ure El we plot the minimum momentum kick Sp — p — k oq 
required to reach the unstable region as a function of 
SI. Here, fc eq is the equilibrium value of the momentum 
ki (fc cq = koV 1 - (mn/2k%) 2 for SI < 4.0fcg/2m and 
fc oq = for S7 > 4.0k 2 /2m). The blue solid and green 
dashed lines correspond to the value of Sp required to 
reach the energetically unstable region for positive and 
negative kicks, respectively. The red dotted line in Fig. [3] 
shows the critical momentum kick needed to reach the 
dynamically unstable region. 

Let us finally notice that, as anticipated in the intro- 
duction, the critical velocity of the supercurrent differs 
from the value predicted using the Landau criterion be- 
cause of the breaking of Galilean invariance [3] . Differ- 
ently from the supercurrent instabilities that can exhibit 
either energetic and dynamical instablty, the Landau in- 
stability has always an energetic nature. 

The dipole oscillation. — We expect that the instability 
discussed above can be observed experimentally through 
the excitation of the collective dipole oscillation in a har- 
monic trap. Actually in our systems the center of mass 
motion is coupled to the spin degrees of freedom by the 
Raman coupling and Kohn's theorem is no longer appli- 
cable (2l|. The dipole mode corresponds to an oscilla- 
tion of the condensate around one of the minima of the 
single-particle spectrum (Fig.[IJa)) and during the oscil- 
lation the condensate behaves as the supercurrent state 




FIG. 3: Critical momentum kick Sp to the condensate which 
causes instabilities, in units of fco as a function of Q, in units 
of kg /2m. The blue solid line and the green dashed line show 
the critical momentum to reach the energetic instability for 
the lowest energy state and the metastable state, respectively, 
when gn — 0.48fcg /2m. The red dotted line shows the critical 
momentum kick required to get the dynamical instability. 



discussed above with alternating velocity. As one can see 
from Fig. [3l if we provide a sufficiently large momentum 
kick to the condensate the system will enter a dynami- 
cally unstable region and the oscillation will break down. 
Indeed in [l7j . where the dipole oscillation was excited 
by suddenly transferring a momentum kick to the con- 
densate, the authors could not observe a single-frequency 
oscillation in the region 2.5/cq/2to < SI < 4.0fco/2m. Our 
estimate from Fig. [3] shows that if the initial momentum 
kick is \Sp\ ps 0.4/c , the condensate actually becomes dy- 
namically unstable in the above region. 

Observing the energetic instability is, in general, more 
challenging than observing the dynamical instability be- 
cause thermalization and collisional effects are required 
to trigger it. Furthermore the corresponding critical ve- 
locity depends on the value of the density and becomes 
smaller and smaller as one approaches the border of the 
condensate. 

In conclusion we have shown that a remarkable conse- 
quence of spin-orbit coupling is the occurrence of dynam- 
ical instability in the superfluid flow when its velocity 
reaches a critical value. Contrary to usual superfluids, 
this unexpected feature takes place also for configura- 
tions with uniform density and is the consequence of the 
violation of Galilean invariance caused by the spin-orbit 
coupling. 
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